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Abstract 

We derive a version of the virial theorem that is applicable to diatomic planetary 
atmospheres that are in approximate thermal equilibrium at moderate temperatures 
and pressures and are sufficiently thin such that the gravitational acceleration can be 
considered constant. We contrast a pedagogically inclined theoretical presentation with 
the actual measured properties of air. 

In his widely discussed article, Miskolczi [T] postulates that the virial theorem, which 
relates the average kinetic and average potential energies of a bound mechanical system (see 
[2] for a thorough introduction), can be applied to a planetary atmosphere in equilibrium in 
the planet's gravitational field. 

To investigate if Miskolczi's postulate is correct (whether or not the postulate was cor- 
rectly applied in pQ is a question beyond the scope of the present paper), let us first consider 
the case of a bouncing ball in an homogeneous gravitational field. (The gravitational field 
within the Earth's atmosphere, the thickness of which is small compared to the Earth's ra- 
dius, is approximately homogeneous. However, the same results presented here could also 
be obtained using a Newtonian gravitational potential [20.) 

Consider dropping the ball from a height hb above the surface, and assume that it falls 
without air resistance, and bounces back from the ground with no loss of kinetic energy. We 
wish to calculate its average kinetic and average potential energy. It is sufficient to compute 
these averages for the first part of the ball's motion, as it falls to the ground; the bounce-back 
is just a time-reversed copy of its initial drop, and afterwards, in the absence of dissipative 
losses, the ball repeats the same motion ad infinitum. 

We presume that the ball was dropped at t — 0. At any other time t > before the ball 
hits the ground, its height will be 

h = h b - l -gt 2 , (1) 
where g (~ 9.81 m/s 2 on the Earth) is the surface gravitational acceleration. From this, 



1 The author wishes to thank Gyorgy Major for bringing Ref. [3] to his attention. 
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and notably, the time it takes to reach the ground (h = 0) is 



^=\l 2 ~J- (3) 

The velocity of the ball at time t (0 < t < to) is 

v = gt= y/2g{h b - h). (4) 
The kinetic energy K and potential energy U of the ball are calculated the usual way: 

ii =-mu = -™? t , (5) 
[/ =mgh = mgh b — -mg 2 t 2 . (6) 
According to the virial theorem, for a bound mechanical system with kinetic energy K, 

2<K) = (Vr a .F a Y (7) 

where () denotes time averaging and r a and F a are the position of, and the force acting on, 
the a-th particle that constitutes the system. In the case of a system in which the potential 
energy is an homogeneous function of degree k of the coordinates, we get 

2(K) = k(U). (8) 

For (J6j), /c = 1. The time averages of the kinetic energy (jSJ) and potential energy ([6]) between 
t — and t — to can be calculated as 

(K) =i / ° K = \mghl = \mgh b , (9) 
*o Jo 6 3 

1 f to 1 2 

(?7) =— / U = mgh b - -mgh b = -mgh b . (10) 
to Jo 

Hence, 2{K) = (U) and the virial theorem for the potential (j6]) appears satisfied. This was 
made possible, in part, by referencing the height h to ground level; this allowed us to ignore 
the effects of the ground surface in ([7]), because at the ground, r a = 0, as observed also in 

®- 

How about a column of atmospheric gas? We assume a column of gas in hydrostatic 
equilibrium standing over a unit surface area in an homogeneous gravitational field. Its 
density will be a function of height h above ground: 

p = p(h). (11) 

The pressure at h is equal to the weight of gas situated at heights above h: 

POO 

p(h) = / gp(h') dti. (12) 
Jh 



We assume that the gas is in thermal equilibrium (the real atmosphere isn't, but that's 
another story), so its temperature is constant: 



T = T . (13) 

We also assume that the gas obeys the ideal gas law (this is a valid approximation for air at 
room temperature and sea level pressure), hence 

pV = nRT, (14) 

where V is the volume of n moles of gas, and R ~ 8.31 JK _1 mol _1 is the ideal gas constant. 
The mass of n moles of gas is nM n where M n (~ 0.029 kg/mol for air) is the molar mass of 
the gas; its density is p = nM n /V, hence V = nM n /p. We can thus rewrite the ideal gas 
law in the form 



M n 



Using this in f lT2|) . we obtain 

RT 

M n 

or, in differential form, 
which can be solved trivially: 



h 



p = I gp(h') dh' - I gp(h') dh', (16) 
RTq dp 



gM n dh 
RTp 1 

gM nP 

RT 



-P, (17) 



dp = dh, (U 



\ogp = h + C, (19) 



or 



gM n 

p = poe -9M n h/RT ^ (2Q) 



which we can also write in the form 



p=p Q e- h ' h \ (21) 
* m1, (22) 



gM n 



For air at T = 273 K, we get 



8.31 x 273 

ho = ~ 8 km, (23) 

9.81 x 0.029 v ; 



which agrees well with the observed properties of the atmosphere. 

So what about the virial theorem? Going back to the bouncing ball for a moment, we 
can immediately spot a potential problem: what if the ball is moving horizontally as well? 
Indeed, it can move horizontally at an arbitrary velocity, yet its potential energy will be 
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no different, hence the virial theorem fails. We must make sure that we only consider the 
vertical component of the ball's velocity before the theorem can be considered valid. The 
velocity of the ball can be written in rectilinear form as v 2 = v\ + v 2 + v 2 , but we are only 
interested in the vertical component. In the specific case when the average velocities in the 
x, y, and z direction are the same, we get {v 2 z ) = |(f 2 )- Accordingly, the virial theorem in 
this case reads 

\{K) = (U). (24) 

This result can also be obtained using another argument, presented in [3]: rather than 
allowing the coordinates to remain unbounded in the horizontal plane, we can consider 
confining the ball to within a cylinder of unit radius, integrating and time averaging the 
force acting on the ball as it hits the cylinder walls, in order to obtain the right-hand side 
of©. 

An atmosphere, unfortunately, is not made of bouncing balls, however appealing that 
picture might appear. Air, in particular, is composed mainly of diatomic gases (notably 
N 2 and O2), which at room temperature have two rotational degrees of freedom in addition 
to the three translational degrees of freedom discussed above. (At higher temperatures, 
vibrational modes also play a role.) The kinetic energy of a column of gas is its internal 
thermal energy. The principle of equipartition of energy states that internal thermal energy 
is distributed equally between all degrees of freedom. Therefore, the virial theorem now 
reads 

\{K) = (U). (25) 

This is our main result, valid for any diatomic atmosphere that obeys the ideal gas law in 
an homogeneous gravitational field at moderate temperatures. 

For a column of gas over a unit surface area, the thermal kinetic energy is 

poo poo 

K= c v T p dh = c v poT / e' h/ho dh = c v p T h , (26) 
Jo Jo 

where cy is the specific heat of the atmosphere at constant volume. The potential energy is 
just the gravitational potential energy: 

poo poo 

U= gphdh = gp / he~ h/ho dh = gp h 2 . (27) 
Jo Jo 

The ratio of the two is 

£ = 4" (28) 
K c v T 

Using f[2~2l . we obtain 

U - = (29) 
K c v M n 1 ' 

Given ( 125|) and ( 1290 . we can calculate the specific heat cy. For air, we obtain 

5R J 
c v = — — ~ 716 — — . 30 
2M n K-kg y 1 
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a value that agrees well with the known properties of air (cy = 718 JK _1 kg _1 ). 

In this derivation, we assumed that T = T is constant. Our result, however, remains valid 
even when T is not constant, so long as the gas is in "local thermodynamic equilibrium", 
which ensures that the principle of equipartition remains valid and that thermodynamic 
quantities, such as the temperature or specific heat, remain well-defined. To see this, we 
first rewrite the condition of hydrostatic equilibrium ffl2l) in differential form: 

dp = —gp dh. (31) 

This allows us to write the potential energy of the gas (1271) as 



OO POO POO J m POO 



U = I gph dh = — I h dp = — I h-^- dh = I p dh, (32) 



o Jo Jo dh j Q 

where the last step was taken by integrating in parts and using p(oo) = 0. On the other 
hand, the thermal kinetic energy (126]) can be rewritten as 

K = c v I" Tp dh = Hp dh. (33) 

Jo R Jo 

The ratio of (132]) and ( 133]) remains the same constant ratio (129]) that we obtained in the 
T = T case: 

U R 

K = ^M„- (34 > 

even as T varies with altitude. Therefore, even as we allow T to be a function of h, (|25j) 
remains satisfied. 

Hence we were able to demonstrate without having to invoke concepts such as "hard 
core" potentials or intramolecular forces that the virial theorem is indeed applicable to 
the case of an atmosphere in hydrostatic equilibrium. However, it must be "handled with 
care": the nature of the atmosphere and the fact that the horizontal (translational) and 
internal (rotational) degrees of freedom of the gas molecules are unrelated to the gravitational 
potential cannot be ignored. 
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